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We report the observation of plasma oscillations in an ul- 
tracold neutral plasma. With this collective mode we probe 
the electron density distribution and study the expansion of 
the plasma as a function of time. For classical plasma con- 
ditions, i.e. weak Coulomb coupling, the expansion is domi- 
nated by the pressure of the electron gas and is described by a 
hydrodynamic model. Discrepancies between the model and 
observations at low temperature and high density may be due 
to strong coupling of the electrons. 



One of the most interesting features of neutral plasmas 
is the rich assortment of collective modes that they sup- 
port. The most common of these is the plasma oscillation 
JIJ , in which electrons oscillate around their equilibrium 
positions and ions are essentially stationary. This mode 
is a valuable probe of ionized gases because the oscillation 
frequency depends solely on the electron density. 

In an ultracold neutral plasma as reported in Q, the 
density is nonuniform and changing in time. A diag- 
nostic of the density is thus necessary for a variety of 
experiments, such as determination of the three-body re- 
combination rate at ultralow temperature || , and obser- 
vation of the effects of strong Coulomb coupling B in a 
two-component system. A density probe would also aid 
in the study of the evolution of a dense gas of cold Ry- 
dberg atoms to a plasma || , which may be an analog of 
the Mott insulator-conductor phase transition 0] . 

In this work we excite plasma oscillations in an ultra- 
cold neutral plasma by applying a radio frequency (rf) 
electric field. The oscillations are used to map the plasma 
density distribution and reveal the particle dynamics and 
energy flow during the expansion of the ionized gas. 

The creation of an ultracold plasma has been described 
in A few million metastable xenon atoms are laser 
cooled to approximately 10 /iK. The peak density is about 
2 x 10 10 cm -3 and the spatial distribution of the cloud is 
Gaussian with an rms radius er ps 220 /im. These parame- 
ters arc determined with resonant laser absorption imag- 
ing (t]]. To produce the plasma, up to 25% of the atoms 
are photoionized in a two-photon excitation. Light for 
this process is provided by a Ti:sapphire laser at 882 nm 
and a pulsed dye laser at 514 nm (10ns pulse length). 
Because of the small electron-ion mass ratio, the result- 
ing electrons have an initial kinetic energy (E e ) approxi- 
mately equal to the difference between the photon energy 
and the ionization potential. In this study we vary E e /ks 



between 1 and 1000 K. The initial kinetic energy of the 
ions varies between 10 /iK and 4mK. 

For detection of charged particles, a small DC field 
(about lmV/cm) directs electrons to a single channel 
electron multiplier and ions to a multichannel plate de- 
tector. The amplitude of the rf field that excites plasma 
oscillations, F, varies between 0.2 — 20mV/cm rms. All 
electric fields are applied to the plasma with grids located 
above and below the laser-atom interaction region. 

In the absence of a magnetic field, the frequency of 
plasma oscillations is given by f e = (l/27r) yj e 2 n e / e^m e 
jij. Here, e is the elementary charge, n e is the electron 
density, eo is the permittivity of vacuum, and m e is the 
electron mass. This relation is most often derived for an 
infinite homogeneous plasma, but it is also valid in our 
inhomogeneous system for modes which are localized in 
regions of near resonant density. Corrections to f e due 
to finite temperature Q depend on the wavelength of 
the collective oscillation, which is difficult to accurately 
estimate. Such corrections are not expected to be large 
and will be neglected. We observe plasma oscillations 
with frequencies from 1 to 250 MHz. This corresponds 
to resonant electron densities, n r , between 1 x 10 4 cm~ 3 
and 8 x 10 s cm -3 . The oscillation frequency is sensitive 
only to n e , but, as explained in the core of the plasma 
is neutral. This implies that plasma oscillations measure 
electron and ion densities in this region (n e = rii = n) . 

Figure [l]a shows electron signals from an ultracold neu- 
tral plasma created by photoionization at time t = 0. 
Some electrons leave the sample and arrive at the de- 
tector at about 1 /is, producing the first peak in the sig- 
nal. The resulting excess positive charge in the plasma 
creates a Coulomb potential well that traps the remain- 
ing electrons In the work reported here, typically 
90 — 99% of the electrons are trapped. Debyc shielding 
maintains local neutrality inside a radius r e beyond which 
the electron density drops to zero on a length scale equal 
to \d- The value of r e depends on the fraction of elec- 
trons that has escaped, and Xo is the Debye screening 
length, \d = \J £okBT e /e 2 n e , where T e is the electron 
temperature. For our conditions r e >, 2<r, and \d <C o. 
As the plasma expands, the depth of the Coulomb well 
decreases, allowing the remaining electrons to leave the 
trap. This produces the broad peak at w 25 [is. 

In the presence of an rf field an additional peak appears 
in the electron signal (Fig. 0a). We understand the gen- 
eration of this peak as follows: The applied rf field excites 
plasma oscillations only where the frequency is resonant. 
Energy is thus pumped into the plasma in the shell with 
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the appropriate electron density [n — n r ). The ampli- 
tude of the collective electron motion is much less than 
a, but the acquired energy is collisionally redistributed 
among all the electrons within 10 — 1000 ns |J, raising 
the electron temperature. This increases the evaporation 
rate of electrons out of the Coulomb well, which produces 
the plasma oscillation response on the electron signal. 

The resonant response at a given time, S(t), is pro- 
portional to the number of electrons in the region where 
the density equals n r . If we make a simple local den- 
sity approximation and neglect decoherence of the oscil- 
lations, S(t) oc F 2 J d 3 r n(r, t) S[n(r,t) - n r \. The width 
in time of the observed signal (Fig. Qa) reflects the den- 
sity distribution of the sample |lfj| . At early times when 
the density is higher than n r almost everywhere, S{t) is 
negligibly small. As the cloud expands and the density 
decreases, the response grows because the fraction of the 
plasma which is in resonance increases. The peak of the 
response appears approximately when the average den- 
sity, fi, becomes resonant with the rf field. S(t) vanishes 
when the peak density is less than n r . 

The resonant response arrives later for lower frequency 
(Fig. |l|b) as expected because n decreases in time. As- 
suming that the plasma density profile remains Gaussian 
during the expansion, S(t) can be evaluated and its am- 
plitude scales as F 2 /n r . In Fig. |l|b the data have been 
normalized by this factor and the resulting amplitudes 
are similar for all conditions. 

By equating h to n r when the response peak arrives, 
we can plot the average plasma density as a function of 
time (Fig. |J). The data are well described by a self similar 
expansion of a Gaussian cloud, n = N/[4iT(a 2 + v 2 t 2 )] 3 ^ 2 , 
where <7o is the initial rms radius and vq is the rms radial 
velocity at long times. N is determined independently by 
counting the number of neutral atoms with and without 
photoionization. The extracted values of cto ar e equal to 
the size of the initial atom cloud. In such an expansion, 
the average kinetic energy per particle is 2>mv 2 /2. 

Figure |3| shows the dependence of vo on density and 
initial electron energy. We first discuss data with E e > 
70 K, for which the expansion velocities approximately 
follow vo = y E e /am,i, where rrii is the ion mass and 
a = 1.7 is a fit parameter. For the plasma to expand 
at this rate, the ions must acquire, on average, a veloc- 
ity characteristic of the electron energy. This is much 
greater than the initial ion thermal velocity. Electron- 
ion equipartition of energy would yield vq = \J E e /3mi, 
close to the observed value. However, due to the large 
electron-ion mass difference, this thermalization requires 
milliseconds Jo) . The observed expansion, in contrast, oc- 
curs on a time scale of tens of microseconds. One might 
expect the expansion to be dominated by the Coulomb 
energy arising from the slight charge imbalance of the 
plasma, but this energy is about an order of magnitude 
less than the observed expansion energy. Also, by Gauss' 
law, it would only be important in the expansion of the 



non-neutral outer shell of the plasma. The oscillation 
probe provides information only on the neutral core be- 
cause it relies on the presence of electrons. 

A hydrodynamic model [jy], which describes the 
plasma on length scales larger than Ap, shows that the 
expansion is driven by the pressure of the electron gas. 
The pressure is exerted on the ions by outward-moving 
electrons that are stopped and accelerated inward in the 
trap. For the hydrodynamic calculation, ions and elec- 
trons are treated as fluids with local densities n a (r) and 
average velocities u a (r) = (v a (r)). Here, a refers to ei- 
ther electrons or ions, and (• ■ ■} denotes a local ensemble 
average. Particle and momentum conservation lead to 
the momentum balance equations 
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Here n a kBT a represents a scalar pressure p!l| . The ion 
and electron equations are coupled by R a h, which is the 
rate of momentum exchange between species a and b. 
The exact form of this term is unimportant for this study, 
but R a f, = — Rha- Plasma hydrodynamic equations typi- 
cally have electric and magnetic field terms, but applied 
and internally generated fields are negligible when de- 
scribing the expansion. 

We can make a few simplifying approximations that 
are valid before the system has significantly expanded. 
The directed motion is negligible, so we set u a ss ev- 
erywhere. Because n e w rii — n, du e /dt « dui/dt. Due 
to the small electron mass, the rate of increase of aver- 
age electron momentum is negligible compared to that 
of the ions. The electron momentum balance equation 
then yields V(nfcsT e ) sa R ei , which describes a balance 
between the pressure of the electron gas and collisional 
interactions. This is the hydrodynamic depiction of the 
trapping of electrons by the ions. 

In the ion momentum balance equation, we eliminate 
Ri e using the electron equation, and we drop the pres- 
sure term because the ion thermal motion is negligible. 
Thus mindui/dt « — V(nfcBT e ), which shows that the 
pressure of the electron gas drives the expansion Jl2| . 
This result implies that the ions acquire a velocity of or- 
der y^fceT^/mi, which is in qualitative agreement with 
the high E e data of Fig. [| To calculate the expansion 
velocity more quantitatively, one must consider that as 
electrons move in the expanding trap, they perform work 
on the ions and cool adiabatically. The thermodynamics 
of this process is beyond the scope of this study. 

The data in Fig. || indicate that about 90 % of the ini- 
tial kinetic energy of the electrons is transferred to the 
ions' kinetic energy, 3mi«Q/2 = 3E e /2a. This does not 
imply that the temperature of the ions becomes compara- 
ble to E e /kB in this process. For the ions, Uj increases, 
but rrii (| Vi — u,;| 2 }, which measures random thermal mo- 
tion and thus temperature, is expected to remain small. 
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This follows from slow ion-electron thermalization |B[ and 
correlation between position and velocity during the ex- 
pansion 

We now turn our attention to systems with E e < 70 K 
(Fig. ||). They expand faster than expected from an ex- 
trapolation of vq = y^E e /arrii, and thus do not even 
qualitatively follow the hydrodynamic model. A relative 
measure of the deviation is (miV 2 — E e /a)/(E e /a). Fig- 
ure ^ shows that the relative deviation increases with in- 
creasing electron Coulomb coupling parameter [Q, r e = 
(e 2 /47T£ a)/k B T e . Here, a = (47m/3)~ 1/3 is the Wigner- 
Seitz radius, n is the peak density at t — 0, and the 
temperature is calculated by 3kgT e /2 — E e . 

The fact that the relative deviation depends only on 
T e , and that it becomes significant as T e approaches 1, 
suggests that we are observing the effects of strong cou- 
pling of the electrons |l5f] . The hydrodynamic model of 
the plasma is only valid when T e <C 1. When r e >, 1, elec- 
tron and ion spatial distributions show short range corre- 
lated fluctuations that are not accounted for in a smooth 
fluid description GS. Correlations between the ion and 
electron positions would provide the excess kinetic en- 
ergy observed in the expansion by lowering the potential 
energy of the plasma. This satisfies overall energy con- 
servation and it may also explain the systematically poor 
fits of the data for high T e (See Fig.|) . 

Strong coupling is also predicted to alter the relation 
for the frequency f e with which we extract the 

plasma density, size, and expansion velocity. The trend 
of this effect agrees qualitatively with the observed devi- 
ation, but knowledge of the wavelength of the collective 
oscillation is needed for a quantitative comparison. 

Other possible explanations for the deviation are re- 
lated to how the ultracold plasma is created. The 
10 ns duration of the photoionization pulse is long com- 
pared to the time required for electrons to move an in- 
terparticle spacing. Photoionization late in the pulse 
thus occurs in the presence of free charges, which will 
depress the atomic ionization threshold by AEjp rs 
ifc B T e {[(3r e ) 3 / 2 + l] 2 / 3 - 1} || . This effect might in- 
crease the electron kinetic energy by AEjp above what 
has been assumed. However, as shown in Fig. || the cal- 
culated AEip is about an order of magnitude smaller 
than the observed effect. The random potential energy 
of charged particles when they are created may also yield 
a greater electron energy than E e M]. 

High r e (high density and low temperature) conditions 
are desirable for studying the three-body recombination 
rate in an ultracold plasma. The theory po| for this pro- 
cess was developed for high temperature, and is expected 
to break down in the ultracold regime ||]. Measuring or 
setting an upper limit for the recombination rate is not 
possible until the dynamics of high r e systems is un- 
derstood. We are currently studying this problem with 
molecular dynamics calculations. 

We have shown that plasma oscillations are a valuable 



probe of the ion and electron density in an ultracold neu- 
tral plasma. This tool will facilitate future experimental 
studies of this novel system, such as the search for other 
collective modes in the plasma and further investigation 
of the effects of correlations due to strong coupling. 
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FIG. 1. Electron signals from ultracold plasmas created by 
photoionization at t — 0. (a) 3 x 10 4 atoms are photoionized 
and E e /kB = 540 K. Signals with and without rf field are 
shown. The rf field is applied continuously, (b) 8 x 10 4 atoms 
are photoionized and E e /kB = 26 K. For each trace, the rf 
frequency in MHz is indicated, and the nonresonant response 
has been subtracted. The signals have been offset for clarity 
and have been normalized by F 2 /n r . The resonant response 
arrives later for lower frequency, reflecting expansion of the 
plasma. For 40 MHz, n r = 2.0 x 10 7 cm" 3 , and for 5 MHz, 
n r = 3.1 x 10 5 cm" 3 . 
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FIG. 2. Expansion of the plasma for TV = 5 x 10 5 
photoionized atoms. The expansion is well described by 
n — N/[Ait(<Tq + v^t 2 )] 3 / 2 . Horizontal error bars arise from 
uncertainty in peak arrival times in data such as Fig. jjjb. Un- 
certainty in TV is negligible in this data set, but is significant 
for smaller TV. The fits are consistently poor at low E e , as in 
the 3.9 K data. 
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FIG. 3. Expansion velocities, vo, found from fits to data 
such as in Fig. [| The initial average density, no, va ries from 
6 x 10 6 to 2.5 x 10 9 cm -3 . The solid line, v = y/E e /am,i, 
with a = 1.7, is a fit to data with E e /k B > 70 K. The be- 
havior of low E e data is discussed in the text. Uncertainty 
in Vo is typically equal to the size of the symbols. There is a 
0.5 K uncertainty in E e /kB reflecting uncertainty in the dye 
laser wavelength. Note that for E e /kB < 70 K, vo shows a 
systematic dependence on no. 




FIG. 4. Excess expansion energy, AE — rriiV 2 — E e /a, rel- 
ative to E e /a, as a function of F e , the Coulomb coupling 
parameter for the electrons at t = 0. The solid line results 
from equating AE to the predicted suppression of the atomic 
ionization potential in the plasma. Horizontal error bars arise 
from uncertainty in E e . Vertical error bars reflect uncertainty 
in both E e and vq. 
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